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The findings of X-ray and neutron scattering experiments on amorphous systems are interpreted 
within the framework of the theory of EucUdean random matrices. This aUows to take into account 
the topological nature of the disorder, a key ingredient which strongly affects the vibrational spectra 
of those systems. We present a resummation scheme for a perturbative expansion in the inverse 
particle density, allowing an accurate analytical computation of the dynamical structure factor 
within the range of densities encountered in real systems. 

I. INTRODUCTION 

Propagating density fluctuation of macroscopic size (hydrodynamic limit) are known to exist both in ordered and 
in disordered materials. Whereas in the ordered ones those excitations (phonons) persist up to momenta of about 
one tenth of the Debye momentum, the fate of excitations of microscopical size in disordered systems is still a quite 
puzzling issue, both from the theoretical and the experimental point of view. 

Recently, a lot of experimental attention (Buchenau et al. 1986, Benassi et al. 1996, Foret et al. 1996, Masciovecchio 
et al. 1996, Masciovecchio et al. 1997, Masciovecchio et al. 1998, Monaco et al. 1998, Sette et al. 1998, Fioretto et al. 
1999, Ruocco et al. 1999, Sokolov et al. 1999) has been paid to the high frequency dynamics of disordered systems like 
glasses. As a matter of fact, high-resolution inelastic X-ray scattering (IXS) and neutron scattering techniques have 
made accessible to experiment the region where the exchanged external momentum p is comparable to p^, namely the 
momentum where the static structure factor has its first maximum: 

O.lpo<p<po (1) 

The typical momentum p^ can be thought as the generalization of the Debye momentum to the case of disordered 
systems. With these data available for the study of excitations of microscopic size, one is naturally led to ask which 
of the well-known features of phonons survive even in disordered systems. A number of facts have emerged from the 
experiments: 

1. The dynamic structure factor S{p,u}) has a Brillouin-like peak for momenta up to p/po ^ 0.5. This inelastic 
peak is due to the interaction of the external photon (or neutron) with some excitation of the system. A very 
controversial issue is the propagating nature of the excitations (see, for example, Caponi et al. and Foret et al. 
in this same volume) within that range of momenta. As a guideline, one could say that as long as the dispersion 
relation between the position of the peak (in the frequency domain) and the external momentum ujp{p) is nearly 
linear, the excitations are likely propagating. Of course, the study of Anderson localization (see, for instance, 
Parisi 1994) shows that, strictly speaking, this condition is neither necessary nor sufficient, so care should be 
taken. 

2. The peak has a width F, whose dependence on the momentum p within the range (^ has been described by 
means of the scaling law (for a large variety of materials): 

F (X (2) 

Let us remark that, while there seems to be a quite general agreement about the fact that F is not affected by 
changes in the temperature, there is still some debate about the actual value of the exponent a. More precisely, 
at low enough momenta (e.g. light scattering) a is undoubtedly ~ 2, while at higher momenta (e.g. X-ray 
scattering) neither the value a ~ 2, nor the value a ^ 4, describing instead a Rayleigh scattering regime, can 
be ruled out on experimental grounds. 

3. The density of states (DOS) exhibits an excess respect to the Debye behavior {g{uj) cx w^), known as the Boson 
peak. This feature is particularly remarkable for strong glasses. Generally speaking, the typical frequencies of 
the Boson peak are in the region where the relation between the frequency of the Brillouin-peak LUp and the 
momentum p is still linear. 
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4. A secondary peak at frequencies smaller than the Brillouin one develops for large momenta, becoming dominant 
for p/po ~ 0.5. 

A number of basic insights on the spectral properties of glasses have been obtained by means of molecular dynamics 
simulations (Mazzacurati et al. 1996, Dell'Anna et al. 1998, Horbach et al. 1998, Ribeiro et al. 1998, Sampoli et 
al. 1998, AUend et al. 1999, Feldman et al. 1999, Taraskin and Elliot 1999) on systems such as argon, silica and 
water. Among the others, let us point out the fact that, in the glass phase, the high frequency dynamics is very well 
described in the framework of the harmonic approximation. Both dynamic quantities like the dynamic structure factor 
(Ruocco et al. 2000) and thermodynamic observables as the specific heat (Horbach at al 1999) have been shown to 
be correctly described by models where only vibrations around quenched positions are taken into account. In figure 
(|l|) the numerical results obtained by Ruocco et al. (2000) are shown, which illustrate the above-mentioned spectral 
features. 

On the theoretical side, the study of the short time properties of glasses has been attempted from two different 
(and somewhat complementary) points of view. 

• On one hand, it has been shown (Gotze and Mayr 2000) that the Mode Coupling Theory (MCT), which usually 
describes the long time limit of the time correlators, can be modified in the glassy phase so as to describe only the 
excitations around the quenched structure. Within such an approach, S{p^ lo) was computed for a hard-sphere 
glass by means of a generalized hydrodynamics approximation. In this framework, a Brillouin peak with a linear 
dispersion relation [propagating excitations) is obtained up to p/po ^ 0.5, as well as a secondary peak for higher 
values. Interestingly, the peak width T was shown to follow the simple scaling law: 



only in the very low momenta regime p/po < 0.1, whereas a different law applies to the range of momenta 
comparable with the experimental and simulation values (see the figure (8) in Gotze and Mayr 2000) which 
seems hard to reduce to the simple form (^). Moreover, it has been shown (Ruocco et al. 2000) that, resorting 
to the harmonic approximation, the memory function can be simply obtained by means of the eigenvalues and 
eigenvectors of the Hessian matrix. 

• A second approach relies on the study of the statistical properties of random matrices (Metha 1991), since 
within the harmonic approximation the whole dynamical features are encoded in the system's Hessian matrix. 
A crucial point is to distinguish two different classes of systems: 

1. vibrational systems whose disorder has a topological origin, as in glasses (Elliot 1983) 

2. vibrational systems on a lattice, where random couplings constants (Schirmacher et al. 1998, Montagna et 
al. 1999, Kantelhardt et al. 2000, Martin-Mayor et al. 2000, Taraskin et al. 2001) are introduced in order 
to mimic the behavior of real glasses. 

The latter class has recently been studied in some detail (Schirmacher et al. 1998, Taraskin et al. 2001) in order 
to give an insight into the spectral properties, in particular into the Boson Peak. It turns out (Martin-Mayor 
et al. 2000), however, that these models definitely miss the p^ behavior of the peak width. Because of the long 
range order due to the underlying lattice, the general beahviour is given instead by 



even at very low momenta. The above result can be easily obtained within the CPA approximation and can 
be checked by the direct inspection of the eigenvectors. The discrepancy with the low momenta behaviour of 
glasses being due precisely to the lack of topological disorder in lattice models, it seems reasonable to address 
the study of the former class. The problem turns out to reduce to the study of a very special class of random 
matrices, which have been called Euclidean Random Matrices (ERM) (Mezard et al. 1999). Is is worthwhile 
to note that this approach would allow to address, besides the study of the high frequency regime of glasses, 
even the problem of instantaneous normal modes (INM) i.e. the statistical properties of the Hessian matrix of a 
liquid at equilibrium (Wu and Loring 1992, Wan and Stratt, 1994, Keyes 1997, Biroli and Monasson 1999), the 
framework being very similar (Cavagna et al. 1999). Nevertheless, in the following we shall focus only on the 
latter. Recently (Martin-Mayor et al. 2001), the dynamical structure factor S{p,uj) was analytically computed. 
The computation relies on a pertubative expansion, the expansion parameter being the inverse density particle 
1/p. Let us recall what those pertubative results undoubtedly have shown: 



(3) 



TcpA oc p'^ 



(4) 
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— At p QO, the dynamical structure factor is composed of a single delta function, whose position changes 
linearly with the external momentum p at low momenta, representing the undamped propagation of a 
sound wave in an elastic medium. This can be easily understood by considering the infinite number of 
particles per wavelength in that limit. 

— For high but finite densities, the disorder of the position of the particles involved in the propagation causes 
a broadening of the Brillouin peak. In other words, the plane waves are no longer eigenstates of the Hessian, 
hence a finite spreading of the eigenfrequencies involved arises. In principle no analogy with the Rayleigh 
scattering should be expected. As a matter of fact, the perturbative computation, performed up to 
order, shows that at low enough momenta the general behaviour is instead given by: 

Term « (l/p)V- (5) 

Yet at density p = 1 a significant deviation from that behaviour is seen as p/po becomes greater than 0.1, 
in a good qualitative agreement with the results of the MCT approach. 

The perturbative calculation has the problem that both the structure factor and DOS end abruptly at a cut-off 
frequency. Moreover, the experimental densities are of order p 1 or less (in reduced units), hence the ability of 
perturbation theory to catch the features of real systems is at least questionable. 

In order to overcome those difhculties, we present a resummation scheme for the perturbative results (Grigera et al. 
2001) which includes terms at all orders in 1/p. Within this scheme, a very good agreement with numerical results, 
even at large frequencies, is found for values of p > 0.3 (reduced units). 



II. HARMONIC APPROXIMATION 



We shall assume that the particles of our system in the glass phase can only oscillate around their equilibrium 
position, claiming that this is enough to predict the high frequency properties. Here by "equilibrium" we do not 
mean thermodynamic equilibrium (the thermodynamics of glasses is an entirely different issue than the one addressed 
here), but rather mechanical equilibrium (i.e. a position where the forces on all particles are zero). For the sake of 
simplicity, let us consider displacements only along a given direction u: 

Xj{t) ^ x"^ + uifj{t) (6) 

Since the number of equilibrium positions {a^^'^} available to the system is actually infinite, growing exponentially 
with the number of particles N ^ one can ask which one is to take in order to compute the vibrational spectra. The 
answer lies on the so called self averaging hypothesis, which has proved to be correct for many observables in disordered 
systems. Hence, we shall assume that in the thermodynamic limit, the spectra obtained for the different realizations 
of disorder is equal to the one computed by considering the average over the disorder. This statistical approach leads 
to a tremendous simplification of our task, as shown in the following, because the model is defined simply by the 
probability distribution of the random variables {a;^'^}, j — \, . . . , N . 

Within that harmonic framework, the energy of the system is (m is the mass of the particles and O is a frequency 
scale) 

We choose units such that m — 1 and U, — 1. Thus the Hessian matrix 

fc=i 

is an Euclidean random matrix, whose spectral properties we are studying. 

In particular, we are interested both in the dynamic structure factor and in the density of states (DOS). 

• The dynamic structure factor S{p,uj), roughly speaking, gives the spectrum of states "excited" by a plane wave 
with a given momentum p. It can be obtained from the Hessian matrix of the potential (|^) in the classical limit 
and in the one excitation approximation (see for instance Martin-Mayor et al. 2001 for a detailed discussion). 
It reads 
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where \n) are the eigenvectors of the Hessian matrix the eigenfrequencies ujn are the square root of the 
eigenvalues (which are aU positive, see equation(^), T is the temperature, the overhne stands for the average on 



the equihbrium positions {x'p}, while \p) stands for a momentum plane wave {{j\p) = exp[ip • Xj'^]/\/N, where 
|j) represents the vector where the j-th particle has displacement u and the other particles do not move). 

• The DOS describes instead the density of 'all' the vibrational states existing within the system. Interestingly, 
at the level of one-excitation approximation, the DOS can be obtained by: 

ffM^r^^^'Hp-oo,^) (10) 



The above theoretical result has been pointed out only very recently (Martin-Mayor et al. 2001). Although it can 
be checked very easily by means of numerical simulations, we believe it would be also very interesting to check 
it on real systems, using the experimental data at large momenta which arc now available. On the other hand, 
at very high frequencies the one-excitation approximation does not hold, and many-excitations contributions 
should be taken into account, hence the relation ( |l0[) is obviously only an approximation, whose reliability in 
describing real systems is a very interesting matter. Moreover, we shall see in the following that ( p^ ) is crucial 
in obtaining a model independent derivation of the actual small p behaviour of the width of the Brillouin peak, 
holding only for topologically disordered systems. As a matter of fact, the same argument cannot be applied to 
lattice models, because the relation ( p^ ) is not true in topologically ordered models. 

As mentioned, the Hessian is averaged over the disorder, i.e. over all the allowed equilibrium positions. Those 
are clearly distributed in a highly correlated manner, due to the hard-core repulsion and long range attraction of 
the potential. It was shown (Martin-Mayor et al. 2001) that, at the level of the superposition approximation, the 
correlations can be taken into account by doing an uncorrelated average while at the same time suitably renormalizing 
the interaction /: if the k point static correlation function is written as 

g{ri,. . .rk) ^ giri) . . .g{rk), (11) 

then the inverse-density expansion for the dynamic structure factor of a system with correlation function g{r) and 
force F{r) can be obtained from the expansion of a fully uncorrelated system defining a dressed interaction: 

f{r) = g{T)F{T) (12) 

where the value of the "spring" constant F{r) is weighted by the probability to find two centers of oscillations at the 
relative distance r. 



III. NON PERTURBATIVE RESULTS 



Let f{p) be the Fourier transform of f{r). The main 'object' to compute is the resolvent of the Hessian matrix (||), 
which depends on the complex variable z ~ lo"^ + irj. It can be written in the following way: 

G{p,z) = r, (13) 

e{p) = p[/(0) - f{p)] . (14) 

where complex self-energy Ti{j),z) has been introduced. Exploiting the well known relation between resolvent and 
dynamical structure factor (letting aside the prefactor k^Tp^, unessential for the w dependence): 

S''^\p,uj)=^ — 3- limImG(p,tj2 +ir;). (15) 

the following connections between the main features of the dynamical structure factor and the self-energy are estab- 
lished: 
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• The 'bare' dispersion relation e(p), which would give the position of the peak in the elastic medium limit, is 
renormalized by the real part of the self-energy z). This gives u;'^'^""™ (p) , the position of the maximum 
of the structure factor in the frequency domain. Let us note that uj^'^^'-°™(p) is certainly linear for small p, as 
expected. 

• The imaginary part T,"{p, z) computed at the position of the peak w = gives the width, T{p), of the 
S^^\p,cj) by means of: 

E"(p, = a;'-^"°™(p)r(p) (16) 



The self-energy can be obtained as a series in 1/p: the fc'^-order corresponds to k particle-label repetitions when 
calculating the moments of the S^^\p,uj) (the details can be found in Martin-Mayor et al. 2001). It is easy to show 
that the sum of all the cactus diagrams (see figure (||)) is given by the solution of the following integral equation: 

1 r d^Q p(fiQ) - fip - q) , 

pj (27r)^ z - e{q) - T.[q, z) 

Interestingly, Equation(p^) provides a model-independent derivation of the p^ scaling of width of the peak of the 
S'^^\p,Lu). Indeed, the large q contribution to the imaginary part of the integral in equation ( pT[ ) is, because of (|o|) 

d^q 
{2n) 



I]'4(p,z) = -™(A) / j:^U{q) - f{p - q)) . (18) 



where gx{X) is the density of states in the domain of eigenvalues (A = Lu^,gxiu;^) = ^). Ifthespect rum is Debye-like 
we have g\{\) oc A°'^, and it is straightforward to show that ( p^ ) is proportional to uj™^°^"^ {p) p'^ . Then the relation 
( p^ ) implies the scaling: 

r(p)cxp2 (19) 

irrespective of the function /(r). Clearly this is only the large q contribution to the integral, but it has been shown 
(Grigera et al. 2001) that it indeed controls the peak width at small p. 



IV. THE GAUSSIAN CASE 



In the following we shall study the simple case 

/(r)=exphrV(2a2)]. (20) 

This may seem an oversimplification, too distant from any realistic case. However, it is not actually the case, at least 
for small momenta. For the sake of comparison we can see in figure (^) the Fourier transform of the fimction g{r)F{r) 
for Argon at very low temperature (^ lOK) together with that of (|2C). We thus expect our results to be reasonable for 
momenta smaller than the first zero of the Fourier transform of g(r)F{r), which in general is close to the maximum of 
the static structure factor, po. Since the Fourier transform of our force decreases an order of magnitude by po = 2/a, 
we shall take this as our po during the following discussion, and a will be our unit of length. 

We have numerically solved equation (|lj) to find the self-energy for several values of p with the Gaussian choice 
(pO|), thus obtaining the structure factor and the density of states, in the eigenvalue domain. In figure (^) we show 
S^{p, A) for several values of the momentum as obtained from equation (O) for p — 1. We also plot the structure 
factor obtained numerically by the method of moments (Turchi et al. 1982, Benoit 1989, Benoit et al. 1992) (for lower 
momenta, the comparison with the method of moments cannot be done due to finite volume effects, Martin-Mayor et 
al. 2001). A very good agreement with the numerical data is achieved. Note the absence of the secondary peak in 
this model (see below). 

We found that for densities down to p ~ 0.6, the agreement is good for all momenta, in fact good results are 
also obtained for the density of states (see figure (^)). For lower densities, the DOS starts to deviate from a Debye 
behavior at small A, in contradiction with the numerical results. Moreover, below p « 0.31 the approximation gives 
a DOS with an unphysical negative spectrum. The negative spectrum develops continuously as a function of p. Also 
included in figure (|^) is the one- loop result. Notice that even for p = 1, the cactus resummation fails to reproduce 
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the exponential decay of the density of states. This is not unexpected, however, due to the non perturbative nature 
(in Xjp) of this tail (Zee and Affleck, 2000). 

It is interesting to look at the small A limit of g\(X)- As seen in figure (|^) (top) the behavior is very nearly 
Debye (i.e. \/A) for small A. An excess of states relative to the Dcbye case develops for higher A, but in a region of 
eigenvalues well beyond the linear dispersion regime (see bottom of figure (^). Thus this peak, similar to that found 
in (Schirmacher et al. 1998), should not be regarded as a Boson peak. 

Finally, let us look at the scaling of the position and the width of the peak in the frequency domain. In figure (^) we 
show the frequency corresponding to the Brillouin peak as a function of the external momentum p. Let us note that 
a nearly linear dependance, (pointing likely to a propagating excitation) persists up to p/po ~ 0.6 — 0.7. In figure (|^) 
instead we plot r(p), obtained by means of (|l^). As expected, the p^ scaling is found for very small momenta, which 
crosses over to a region where a simple law as (H) is not suitable to describe the real beahviour of the system. 

Note that the region where the scafing is actually found, i.e. p/po < 0.1 is quite different from the region explored 
by X-rays and neutrons scattering experiments, which rather spans the momenta 0.1 < p/po < 0.5, It is worthwile to 
note that the same conclusion can be drawn from the results of Gotze and Mayr using MCT for hard spheres (Gotze 
& Mayr 2000). 

V. CONCLUSIONS AND OUTLOOKS 

In summary, we have applied the random-matrix approach to the study of the high-frequency excitations of glassy 
systems. We have presented a resumation scheme (Grigera et al. 2001) that greatly enhances the predictive power 
of the 1/p expansion prevoisuly obtained (Martin-Mayor et al. 2001). We have compared our analytical calculation 
with numerical results, finding that for not too low densities the only failure of the resumation scheme is its inability 
to repoduce the exponential decay of the density of states. 

Our calculations have been performed choosing a gaussian force, in order to modelize the interparticle force dressed 
with the pair correlation function, and only coUinear displacements have been considered. We believe that our results 
are relevant for realistic glasses, at least in the regime p/po < 1. The basic equation ( [l7| ) shows that, in the general 
case, the position of the spectral peak is linear at small momenta, and the numerical solution of the integral equation 
in the Gaussian case shows that the linearity persists up to p/po ~ 0.6 — 0.7. 

Moreover, the width of the spectral peak turns out to be proportional to p^ in the limit of p ^ (see equation (p^)) 
irrespective of the potential function ,f{x). Interestingly, in the Gaussian case the scaling law holds for a momentum 
range one order of magnitude smaller than the experimental one. At larger momenta, the law is more complicated 
(see figure (0)). Being this result quite similar to MCT result for an hard-spheres system (Gotze and Mayr 2000), we 
believe that it could have a certain degree of univerality. 

This model lacks an important feature of the experimental spectra, namely the secondary peak of the S{q, lo) (which 
some authors have interpreted as the Boson peak, see Horbach et al. 1998). We believe this to be related to two 
important ingredients missing in our model: the vectorial nature of the vibration, and a detailed consideration of the 
particle correlations. The latter should not affect our results for momenta much smaller than the first maximum of 
the static structure factor, pq. Since the secondary peak appears for p ~ poj we should not really expect to describe 
it: in this momentum range the dispersion relation is no longer monotonic in real glasses. Furthermore, Dell'Anna et 
al. 1998 have suggested that transversal excitations may play a prominent role in the raising of the secondary-peak. 
Our approach can be extended to include transverse displacements, though, and work in this direction is in progress. 

VMM was partly supported by CICyT AEN99-0990, AEN97-1693 and M.E.C. TSG was supported in part by 
CONICET (Argentina). 
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FIG. 1. Dynamical structure factor for argon, obtained numerically by Ruocco et al. 2000 at low temperature ~ IQK. The 
Brillouin-like and the secondary peak appear quite clearly. 
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FIG. 2. Cactus diagrams of the 1/p expansion. The numbers correspond to the particle-label repetitions 
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FIG. 4. Dynamical structure factor in the eigenvalues domain for several values of the momentum at p = 1. 
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FIG. 6. Scaling behavior of the position of the peak, in the frequencies domain, for rho = 1. 
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